Zigzag Phase Transition in Quantum Wires 
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We study the quantum phase transition of interacting electrons in quantum wires from a one- 
dimensional (ID) linear configuration to a quasi-lD zigzag arrangement using quantum Monte Carlo 
methods. As the density increases from its lowest values, first, the electrons form a linear Wigner 
crystal; then, the symmetry about the axis of the wire is broken as the electrons order in a quasi- ID 
zigzag phase; and, finally, the electrons form a disordered liquid-like phase. We show that the linear 
to zigzag phase transition is not destroyed by the strong quantum fluctuations present in narrow 
wires; it has characteristics which are qualitatively different from the classical transition. 



Interacting one-dimensional (ID) systems have been 
a fruitful field of study in both condensed matter and 
atomic physics [TJ [2|. Experiments on semiconductor 
quantum wires and carbon nanotubes, for instance, have 
yielded a rich set of data on the ID electron gas over 
the past two decades f^- Linear ion traps, on the other 
hand, provide new systems for studying fundamental ID 
physics [4j as well as potential platforms for quantum 
computing [5] and quantum simulation [6]. These ex- 
perimental systems are not truly ID, of course, and the 
transition from ID to higher-dimensional behavior is of 
both practical and theoretical interest. Here we study 
the first stage in such a transition: the change from a ID 
linear arrangement of particles to a quasi- ID zigzag con- 
figuration, and then to a liquid state at higher densities. 

At low densities, electrons confined to ID by a trans- 
verse harmonic potential form a linear Wigner crystal 
[HIlHS], as illustrated in Fig.jlja). As the electron den- 
sity is increased (or the harmonic confinement relaxed), 
the Coulomb repulsion between particles becomes com- 
parable to the confining potential. The linear crystal 
buckles at a critical value of the electron density [TOl |TT] , 
breaking the symmetry about the longitudinal axis and 
forming a zigzag structure, as depicted in Fig.jljb). This 
system has been studied theoretically in both the weakly 
and strongly interacting limits [9j [T2HT4] : the zigzag 
transition is predicted to be an Ising-type quantum phase 
transition in the strongly interacting limit [12], whereas 
at weak coupling, the critical exponents are non-universal 
[121 [14]. Furthermore, the linear and zigzag phases are 
expected to have only one gapless excitation mode, which 
corresponds to longitudinal sliding of the crystal. At 
higher densities, the integrity of the zigzag structure is 
destroyed, and a second mode becomes gapless. Quasi- 
lD lattice structures were noted in numerical calculations 
on wires with weak confinement . Experimentally, 
evidence for a coupled two-row structure has been ob- 
served in the conductance of quantum wires fabricated 
in GaAs/AlGaAs heterostructures [T6[[T7]. 

Atomic systems provide an alternative to electrons for 
studying zigzag physics. Ions in linear traps, for in- 




FIG. 1. (a) Electrons (red spheres) confined by a harmonic 
potential form a linear Wigner crystal, (b) As electron den- 
sity increases, symmetry about the wire axis breaks and the 
electrons form a zigzag structure. 

stance, have been observed to undergo a zigzag transi- 
tion [18, 19 . Though these experiments are understood 
purely classically [11 , recent theoretical work suggests 
that the quantum zigzag transition may be experimen- 
tally accessible in ion trap systems j2QH2?. Dilute dipo- 
lar gases also undergo a related transition, studied in 
recent theoretical and numerical work [23l [24] . 

To connect to experiments in electronic systems, the 
behavior at intermediate interaction strength is key. We 
observe the zigzag transition at intermediate interaction 
strength, characterizing it through both the pair den- 
sity and a correlation function of the zigzag order. By 
studying the long-range zigzag correlations, we demon- 
strate that the quantum phase transition occurs at pa- 
rameters relevant to quantum wire experiments and is 
qualitatively different from the classical transition. 

The quantum wire consists of TV spinless electrons with 
Coulomb interactions confined to the circumference of a 
ring of radius by a harmonic potential with frequency 
cj, as described by the two-dimensional Hamiltonian 

^ i=l ^ i=l i<j<N^^' 

(1) 

We use effective atomic units: the effective mass m*, 
the electric charge e, the dielectric constant e, and fi are 
all set to 1. In GaAs, the effective Bohr radius = 
h^e/m*e'^ is 9.8 nm, and the energy scale — the effective 
Hartree H* = e^/eag — is 11.9 meV. Using a ring is a 
convenient way to impose periodic boundary conditions 
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on the quantum wire. 

Two length scales in our system are of particular in- 
terest: the Wigner-Seitz radius rg = l/2niD (where niD 
is the linear density), and the length scale vq at which 
the Coulomb interaction between neighboring electrons 
becomes comparable to the harmonic confinement, 

ro ^ f^. (2) 

The zigzag transition occurs when the length scales 
and ro become comparable The classical transition 
has been studied for electrons in liquid helium [10] and 
ion trap systems [TT] . For our system, calculation of the 
classical critical point is straightforward, and the width 
between rows scales as V^n^^^^^criticai [9^ . The classical de- 
scription is valid when the effective Bohr radius is much 
smaller than the interparticle separation at the zigzag 
transition; i.e., when ro ^ 1. Note that for large ro 
(small cj, thus wider wires) the electrons are effectively 
more strongly interacting. For smaller values of ro (nar- 
rower wires), quantum fluctuations play an important 
role. We focus on two values of the confinement, u = 0.1 
and 0.6, which correspond to ro = 5.9 and 1.8, respec- 
tively, because they correspond to experimentally mea- 
sured parameters in GaAs/AlGaAs quantum wires [25j . 
At uj = 0.6 in particular, quantum effects play an impor- 
tant role in the transition since ro is close to 1. 

We calculate the ground state properties of our system 
using Quantum Monte Carlo (QMC) techniques [26|l27]. 
In the first step of our QMC calculation. Variational 
Monte Carlo (VMC), we minimize the variational energy 
of a Slater- Jastrow type trial wave function ^t(I^) = 
J(R)L>(R) [28], using methods described in [29] and [30]. 
We consider three qualitatively different types of single- 
particle orbitals to build the Slater determinant D{IV) — 
localized floating gaussians [31], planewaves, and orbitals 
from density functional calculations — and use the type 
yielding the lowest variational energy at a given density. 
After optimizing the variational parameters, we use Dif- 
fusion Monte Carlo (DMC) to project the trial wavefunc- 
tion onto the fixed-node approximation of the ground 
state [26l[27]. The fixed- node DMC wavefunction is the 
lowest-energy state with the nodes of the trial wavefunc- 
tion and its energy is an upper bound on the true 
ground state energy. We use an extrapolated estimator 
(O)qmc = 2(0)dmc — (O)vMC to calculate observables 
O that do not commute with the Hamiltonian [26 . 

The pair density, defined as (p(r, 6>)p(r^, 0)) /(p(r^,0)) 
where p{r,0) is the local density operator and (r^,0) is 
the location of a fixed electron, is a key microscopic quan- 
tity that allows direct visualization of the system's quasi- 
ID nature. Fig. [2] shows the pair densities in different 
phases of a quantum wire with cj = 0.1. At r^ = 4.0 — the 
"linear phase" — modulations in the pair density indicate 
that the electrons are quasi-localized in a linear arrange- 
ment, as observed in previous QMC calculations [8] [32]. 
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FIG. 2. Pair density, (p(r, 6')p(r^, 0))/(p(r^, 0)), for electrons 
at progressively higher densities in a wire with uj — 0.1. The 
red stars mark r^. (r is plotted relative to f in units of the 
effective Bohr radius , and is plotted in units of the inter- 
particle spacing 2ti/N .) (a) r^ = 4.0, = 30. At low densi- 
ties, the electrons form a linear Wigner crystal, (b) Vs — 3.6, 
= 30. As the density increases, a zigzag structure forms. 

(c) Ts — 3.0, N — 30. The amplitude of the zigzag increases. 

(d) rs = 2.0, = 60. At higher densities, the zigzag structure 
is destroyed. The color scale shows the pair density relative 
to the maximum value for each plot. Only a portion of the 
full periodic system is shown. 



Fig.[2|a) shows that this linear ordering persists even 
when r^ deviates significantly from the center of the wire: 
there is, of course, some short-range zigzag correlation 
in response to the off-axis electron, but at larger dis- 
tances the arrangement is linear. We observe this lin- 
ear phase until the density increases past r^ = 3.79, 
where the zigzag transition occurs; this value is quite 
close to that for the classical transition, rf^^^- = 3.75 for 
Lj = 0.1. Fig.[2|b) shows the system in the "zigzag phase" 
at Tg = 3.6: the electrons are arranged in a long-range 
zigzag pattern. Beyond the transition, the amplitude 
of the zigzag structure continues to increase, reaching 
a maximum value near r^ = 3.0, shown in Fig.[2|c). 

We now turn to the case with stronger quantum fluc- 
tuations, UJ = 0.6, where ro is on the same scale as the 
effective Bohr radius (ro = 1.8). At r^ = 1.5, Fig.jsja) 
shows the system in the linear phase, but modulations in 
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FIG. 3. Pair densities for electrons in a wire with uj — 0.6. 
The red stars mark r^. (r is plotted relative to f, and is 
plotted in units of the interparticle spacing 2n/N .) {3) Vs — 
1.5, A/" = 30. At low densities, electrons in a quantum wire 
form a linear Wigner crystal, (b) Vs = 1.3, = 30. As the 
density increases, a zigzag structure forms, though quantum 
fluctuations smear our correlations in the pair density, (c) 
Vs = 0.5, N = 60. At higher densities, the zigzag structure 
vanishes. The color scale shows the pair density relative to 
the maximum value for each plot. Only a portion of the full 
periodic system is shown. 



the pair density are much smaller than at uo = 0.1. The 
zigzag transition occurs between = 1.4 and Vg = 1.45; 
this deviates significantly from the classical transition 
point, rf""^^ = 1.19 for uj = 0.6. Fig.js^b) shows the 
system in the zigzag phase at = 1.3; quantum fluctu- 
ations have smeared out correlations in the pair density. 
We shall now show, however, that there are strong zigzag 
correlations present despite the rather weak features in 
the pair density. 

To characterize the long-range zigzag order more quan- 
titatively, we first number the electrons along the wire 
axis (i.e., by increasing 0). We can then define a zigzag 
correlation function 



(3) 



where y = r—f denotes the transverse coordinate (for our 
ring geometry, the radial position r relative to the mean 
f). This correlation function corresponds to the field 
{—lyyi that orders in the zigzag state [12]; CzzH'f^ ~ j\) 
indicates how strongly the and j*^ electrons are locked 
in a zigzag pattern as a function of the number of inter- 
vening electrons. Note that Czz is similar to a staggered 
spin correlation function for an antiferromagnetic system. 

Fig.|4] shows the zigzag correlation function, normal- 
ized by the mean squared wire width (y^), plotted at 
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FIG. 4. The zigzag correlation function, normalized by the av- 
erage wire width; - JD/iv^) = ((-l')2/i(-l'')y.>/(y'), 

plotted for various values of rs at (a) cj = 0.1 and (b) u = 0.6 
for = 30 electrons. (^ = r — f is measured in units of the 
effective Bohr radius a^.) In the linear phase (dashed lines), 
Czz decays to 0; in the zigzag phase (solid lines), Czz sat- 
urates to a finite value, indicating long-range zigzag order. 
In the liquid phase (dotted lines) Czz again decays. [Note 
that since our system is periodic, Czz{k) = Czz{N — k), thus 
Czz{k) must be flat at A; = N/2.] 



several values of for uj = 0.1 [Fig.|4ja)] and lj = 0.6 
[Fig.gb)]. In the linear phase (r^ = 3.8,4.0 for a; = 0.1; 
rs = 1.5, 2.0 for a; = 0.6), there is no long-range zigzag or- 
der, and Czz decays to within a few inter-particle spac- 
ings. In the zigzag phase (r^ = 3.0, . . . ,3.7 for uj = 0.1 
and = 1.1, . . . , 1.4 for = 0.6), Czz saturates to a fi- 
nite value. At cj = 0.6, long-range zigzag order is present 
even in the absence of strong long-range positional order 
along the axis of the wire. This is possible because the 
zigzag order is not local (tied to the coordinate along the 
wire axis), but rather depends non- locally on the num- 
bering of the electrons along the wire. This demonstrates 
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that the zigzag transition occurs in the quantum regime. 

When we increase the density further, the zigzag cor- 
relation function again decays, indicating that the zigzag 
structure disappears. This is visible in Fig. [4] at higher 
densities (r^ = 1.0,2.0 for uj = 0.1, = 0.5,0.8 for 
uo = 0.6). At these higher densities, the pair density plots 
show little structure, indicating that the positional order 
has been lost, as seen in Fig.[2|d) at cj = 0.1, = 2.0, 
and in Fig.jsj^c) at cj = 0.6, Vs = 0.5. (We plot results 
from a larger system size, N = 60, in these cases to 
lessen the effects of our ring geometry.) Two rows are 
visible in the pair density, but there are no strong zigzag 
correlations. We identify this liquid-like phase with the 
two-gapless-mode phase described in [12]. 

There are a number of unusual features in Czz- First, 
for values of close to the zigzag transition (e.g., = 3.7 
ai uj= 0.1, and = 1.4 at co'= 0.6), Czz decreases sharply 
before saturating, while the decay is more gradual for 
smaller r^. We also see oscillations in Czz in the phase 
where the zigzag structure is destroyed (r^ = 1.0,2.0 at 
uj = 0.1, and = 0.5,0.8 at = 0.6); these oscillations 
appear to be caused by finite-size effects since we observe 
longer wavelengths at larger N. Finally, at the highest 
densities, Czz shows anti-zigzag ordering [33] at \i—j\=2. 

The long-range value of the zigzag correlation yields 
the order parameter of the phase transition, Mzz- We 
estimate M^^ by averaging over the long-range part of 
the correlation. 
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Mzz is related to the amplitude of the zigzag structure — 
it is half the width between the two zigzag rows. Fig.|5] 
shows the order parameter at = 0.1 and 0.6 as a func- 
tion of (scaled in units of ro). Mzz increases sharply at 
the transition from the linear phase to the zigzag phase. 
As with the pair densities and correlation functions, Mzz 
shows the same generic behavior at = 0.6 as at 0.1 — 
behavior consistent with that of a continuous phase tran- 
sition. The transition occurs at a scaled density close to 
the classical value for cj = 0.1 but at a considerably lower 
scaled density for u = 0.6. In both cases, the behavior 
near the transition differs qualitatively from the classical 
case [9 , highlighting its quantum nature. As the system 
evolves at higher density from the zigzag to the liquid-like 
phase, Mzz decreases gradually. 

In conclusion, we have demonstrated that the linear 
to zigzag transition occurs at parameters relevant to ex- 
periments in quantum wires. A transition to a phase 
with long-range zigzag correlations occurs even in narrow 
wires where large quantum fluctuations smear out den- 
sity correlations; the quantum phase transition in these 
wires differs substantially from the classical case. Future 
work will examine other signatures of this transition; in 
addition, the role of spin and finite temperature at inter- 
mediate interaction strength remains unexplored. 
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FIG. 5. The "Zigzag Amplitude" order parameter Mzz, as 
a function of rs at a; = 0.1 (blue circles) and a; = 0.6 (red 
squares). As is decreased beyond the critical value and 
the system enters the zigzag regime, Mzz increases sharply; 
at lower rs, Mzz decreases gradually as the system enters the 
liquid regime. The behavior deviates significantly from the 
y/n — ricriticai bchavior for classical electrons (black dashed 
line). Dotted lines are a guide to the eye. (Lengths are scaled 
by ro . We note that there are points for both values of uj where 
Mzz seems to be lower than one would expect by drawing 
a curve through the other points; we believe that this is a 
systematic error from the VMC optimization step of our QMC 
calculation.) 
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[33 1 This is due to the decoupling of the two rows; at high 
densities, there are likely to be configurations where two 
"adjacent" (ordered by their coordinate) electrons are 
in the same row, and the next electron is in the opposite 
row. (Le., the \i — j\ =0 and \i — j\ = 1 electrons are 
in the same row, and the \i — j\ = 2 electron is in the 
opposite row, or both the \i — j\ = 1 and \i — j\ = 2 
electrons are in the row opposite the — = electron.) 
Such configurations will give a negative contribution to 

a.(N-ji = 2). 



